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ABSTRACT 

A bright quasar residing in a dense and largely neutral intergalactic medium (IGM) at high redshifts (z > 6) 
will be surrounded by a large cosmological Stromgren sphere. The quasar's spectrum will then show a sharp 
increase in resonant Lyman line absorption at wavelengths approaching and shorter than that corresponding to the 
Stromgren sphere's boundary along the line of sight. We show here that simultaneously considering the measured 
absorption in two or more hydrogen Lyman lines can provide the dynamical range required to detect this feature. 
We model broad and robust features of the Lyman a and Lyman (3 regions of the spectrum of the z = 6.28 quasar 
SDSS J1030+0524, using a hydrodynamical simulation. From the steep wavelength-dependence of the inferred 
absorption opacity, we detect the boundary of the Stromgren sphere at a proper distance of 6.0 ± 0.2 Mpc away 
from the source redshift. From the spectrum alone, we also find that beyond this distance, cosmic hydrogen 
turns nearly neutral, with a neutral fraction of xhi > 0.2, and that the ionizing luminosity of this quasar is in the 
range (5.2 ±2.5) x 10 56 photons sec" 1 . The method presented here, when applied to future quasars, can probe the 
complex topology of overlapping ionized regions, and can be used to study the details of the reionization process. 
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Subject headings: cosmology: theory - cosmology: early universe - galaxies: formation - galaxies: 
high-redshift - galaxies: quasars: general - galaxies: quasars: absorption lines 



1. INTRODUCTION 

The ionization state of the IGM at redshift 6 < z < 7 has been 
a subject of intense study over the past few years. While stud- 
ies of the cosmic microwave background anisotropics by the 
Wilkinson Microwave Anisotropy Probe (WMAP) satellite im- 
ply that the IGM is significantly ionized out to z ~ 15, several 
pieces of evidence suggest that it has a large neutral fraction at 
z ~ 6-7 (see, e.g., Haiman 2003 for a recent review). If in- 
deed the intergalactic hydrogen is largely neutral at z ~ 6, then 
quasars at this redshift should be surrounded by large ionized 
(HII) regions (Madau & Rees 2000; Cen & Haiman 2000). The 
imprint on the absorption spectrum of the damping wing of ab- 
sorption by neutral hydrogen outside the HII region (Mesinger 
et al. 2004), and the size of the HII region itself (Wyithe & 
Loeb 2004) can be used to place stringent limits on the neutral 
fraction of the ambient IGM. 

In principle, the quasar's absorption spectrum contains a record 
of the neutral fraction as a function of position along the line of 
sight. In Figure 1 we illustrate a model for the optical depth 
to Lyman a absorption as a function of wavelength towards 
a zq = 6.28 quasar embedded in a neutral medium (xhi =1), 
but surrounded by a Stromgren sphere with a comoving ra- 
dius of Rg = 44 Mpc. Around bright quasars, such as those 
recently discovered (Fan et al. 2001, 2003) at z ~ 6, the proper 
radius of such Stromgren spheres is expected to be R$ ~ 7.7 

*ra /3 (Nq/6.5 x 10 57 s" 1 ) 1 ' 3 (t Q /2 x 10 7 yr) 1/3 [(l+z G )/7.28]" 1 
Mpc (Madau & Rees 2000; Cen & Haiman 2000). Here x m is 
the volume averaged neutral fraction of hydrogen outside the 
Stromgren sphere and Nq, tQ, and zq are the quasar's produc- 
tion rate of ionizing photons, age, and redshift. The fiducial val- 
ues are those estimated for J1030+0524 (Haiman & Cen 2002; 
Wyithe & Loeb 2004). The mock spectrum shown in Figure 1 
was created by computing the Lyman a opacity from a hydro- 
dynamical simulation (the procedure is described in § 2 below). 
The optical depth at a given observed wavelength, A b s , can be 
written as the sum of contributions from inside (t#) and out- 



side (tq) the Stromgren sphere, TL ya = tr + to- The residual 
neutral hydrogen inside the Stromgren sphere at redshift z < Zq 
resonantly attenuates the quasar's flux at wavelengths around 
A a (l +z), where X a = 1215.67 A is the rest-frame wavelength 
of the Lyman a line center. As a result, tr is a fluctuating func- 
tion of wavelength (solid curve), reflecting the density fluctu- 
ations in the surrounding gas. In contrast, the damping wing 
of the absorption, td, is a smooth function (dashed curve), be- 
cause its value is averaged over many density fluctuations. As 
the figure shows, the damping wing of the absorption from the 
neutral universe extends into wavelengths A b s > 8720 A, and 
can add significantly to the total optical depth in this region. 

The sharp rise in tq at wavelengths A b s < 8720 A is a unique 
feature of the boundary of the HII region, and corresponds to 
absorption of photons redshifting into resonance outside of the 
Stromgren sphere. The detection of this feature has been re- 
garded as challenging: since the quasar's flux is attenuated by 
a factor of exp(-TL yQ ), an exceedingly large dynamical range 
is required in the corresponding flux measurements. However, 
we show here that simultaneously considering the measured ab- 
sorption in two or more hydrogen Lyman lines can provide the 
dynamical range required to detect this feature. In particular, 
we model broad features of the Lyman a and Lyman (3 regions 
of the absorption spectrum of the z = 6.28 quasar SDSS J 1030+0524. 
We find robust evidence for the boundary of the Stomgren sphere, 
and derive new limits on the neutral fraction of the ambient 
IGM, and on the ionizing emissivity of the quasar. 

'Throughout this paper we assume a standard cold-dark matter cosmology 
(ACDM), with (n Al n„, Cl b , n, c% H ) = (0.73, 0.27, 0.044, 1, 0.85, 71 km 
s Mpc~' ), consistent with the recent results from WMAP (Spergel et al. 2003). 
Unless stated otherwise, all lengths are quoted in comoving units. 
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FIG . 1 . — Model from a hydrodynamical simulation for the optical depth con- 
tributions from within (77;) and from outside (to) the local ionized region for 
a typical line of sight towards a zq = 6.28 quasar embedded in a fully neutral, 
smooth IGM, with Rs = 44 Mpc and f mrt = 1 . The dashed curve corresponds 
to to, and the solid curve corresponds to tr. The total Lyman a optical depth 
is the sum of these two contributions, TL ya =tr + to- The dashed-dotted lines 
demarcate the three wavelength regions used for our analysis described in the 
text. In our analysis, the optical depths are averaged over 3.5 A wavelength 
bins, which decreases the fluctuation of tr. For reference, the redshifted Ly- 
man a wavelength is at 8852 A, far to the right off the plot. 

2. ANALYSIS 

The observational input to our analysis is the deepest avail- 
able absorption spectrum of SDSS J1030+0524 (White et al. 
2003). The flux detection threshold in the Lyman a and Ly- 
man (3 regions of this spectrum correspond to Lyman a optical 
depths of Tij m ( LyQ ) w 6.3 and Tii m(Ly/ 3) ss 22.8 respectively. 2 The 
ratio of these two numbers is smaller than that of (f a \a)/(fp Xp), 
where / is the oscillator strength, mainly because of the need to 
subtract foreground Lyman a absorption in the Lyman (3 region. 
The Keck ESI spectrum of SDSS J1030+0524 exhibits a strong 
Lyman a Gunn-Peterson (GP) trough, with no detectable flux 
between wavelengths corresponding to redshifts 5.97 < z < 6.20, 
as well as a somewhat narrower Lyman (3 trough between 5 .97 < 
z < 6.18 (White etal. 2003). 

To summarize these constraints, we have divided the spec- 
trum into three regions, shown in Figure 1. In Region 1, with 
A bs > 8752.5 A, the detection of flux corresponds to the upper 
limit on the optical depth TL ya < 6.3. Region 2, extending from 
8725.8 A < A bs < 8752.5 A, is inside the Lyman a trough, but 
outside the Lyman (3 trough. Throughout this region, the data 
requires 6.3 < TL yQ < 22.8. Region 3, with A b s < 8725.8 A 
has a lower limit TL ya > 22.8. As defined, each of these three 
regions contains approximately eight pixels. 

2 For our purposes, these optical depths can be taken as rough estimates. 
Their precise values are difficult to calculate, with Tij ra (L y/ 3) especially uncer- 
tain (Songaila 2004; Lidz et al. 2002; Cen & McDonald 2002; Fan et al. 2002). 
However, we have verified that our conclusions below remain unchanged when 
the threshold opacities are varied well in excess of these uncertainties. In partic- 
ular, considering ranges as wide as 5.5 < Tij m (Ly Q ) < 7 and 10 < Ti; ra (L y/ 3) < 30 
would lead to constraints similar to those we derive below. 



We modeled the absorption spectrum, attempting to match 
these gross observed features. We do not perform a far more 
demanding pixel-by-pixel fit to the data, given the relatively 
poor quality of the spectra currently available. We utilized a 
hydrodynamical simulation that describes the density distribu- 
tion surrounding the source quasar at z = 6.28. The details of 
the simulation are described elsewhere (Cen 2004, in prepara- 
tion; Mesinger et al. 2004). We extracted density and velocity 
information from 100 randomly chosen lines of sight (LOSs) 
through the simulation box. 3 Along each line of sight (LOS), 
we computed the Lyman a absorption as a function of wave- 
length. The size of the ionized region (Rs) and the fraction of 
neutral hydrogen outside it (%) were free parameters. In ad- 
dition, the quasar's ionizing luminosity, L„, was a third free 
parameter, f loa , defined by L v = /i on U v . Here Lf u = 1 .55 x 10 31 
(v/vhT 1 * [(1 +z)/(l +Zq)]-°\ v h = 3.29 x 10 15 Hz is the ion- 
ization frequency of hydrogen, and zq is the source redshift 
(zq = 6.28 for SDSS J1030+0524). U v results from redshift- 
ing a power-law spectrum with a slope of v~ x 8 , normalized 
such that the emission rate of ionizing photons per second is 
1 .3 x 10 57 , matching the emissivity one obtains (Haiman & Cen 
2002) for SDSS J1030+0524 using a standard template spec- 
trum (Elvis et al. 1994). This translates to an emission rate 
of ionizing photons of Nq « /i on x 1.3 x 10 57 s" 1 . To under- 
stand the significance of these parameters, note that changing 
Rs moves the dashed (td) curve in Figure 1 left and right, while 
changing xui moves it up and down; changing fi on moves the 
solid (tr) curve up and down. Values of xhi < 1 imply the exis- 
tence of an ionizing background flux, with an ionization rate of 
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x 10 s . We add this background flux to that 



of the quasar (the latter dropping as r~ 2 with distance), some- 
what reducing the neutral hydrogen fraction inside the ionized 
region. 

We evaluated tr and To for each LOS, in the range of ob- 
served wavelengths 8700 A < A b s < 8770 A, and smoothed 
the optical depth by averaging exp(-TL ya ) over ~ 3.5 A wave- 
length bins, corresponding to the Keck ESI spectral resolution 
of R = 2500 (White et al. 2003). This procedure was repeated 
for each combination of Rs, Xm, and fi on . We allowed xhi to 
range from 1 to 10~ 3 , to trace its currently allowed range (Fan 
et al. 2002). Values of R s < 42 Mpc (5.8 Mpc proper) are im- 
mediately ruled out by the presence of flux at z > 6.18 in the 
redshifted Lyman (3 region of the spectrum. Hence, we let Rs 
range from 42 to 54 Mpc, in increments of 1 Mpc. Finally, /; on , 
was varied from 0.1 to 10, in increments of 0.1. For each point 
in the three-dimensional parameter space of Rs, Xm, and f loa , 
we computed the fraction of the 100 LOSs that were accept- 
able descriptions of the spectrum of SDSS J1030+0524, based 
on the criteria defined above. To be conservative, we allowed 
up to two pixels in each region to lie outside the allowed range 
of Lyman a optical depths; each LOS that had more than two 
pixels fail the above criteria in any of the three regions was re- 
jected. We find that a more stringent criterion of allowing only 
a single 'faulty' pixel would strengthen our conclusions. 

3 The density and velocity distributions are biased near the density peaks 
where a quasar may reside. However, the spectral region of interest lies well 
outside these biased regions, which, in the case of high-redshift quasars, ex- 
tends on average only to < 1 proper Mpc(Barkana & Loeb 2004). The use of 
randomly chosen LOSs is therefore an accurate statistical representation of the 
expected density and velocity fields. 




FIG. 2. — The percentages of LOSs that pass our analysis criteria. Results 
for three different values of the ionizing flux parameter, f mrt = 0.2, 0.4, 0.6 
are shown in the top, middle and bottom panels, respectively. The fraction of 
accepted LOSs peaks sharply at a maximum of 24% at the parameter choice of 
(Rs, xhi, /ion) = (42 Mpc, 1, 0.4). The radius of the ionized region surrounding 
the quasar is limited to 42 Mpc < Rg < 47 Mpc. The ionization flux parameter 
is tightly constrained around f mn = 0.4, with no LOSs passing our criteria with 
values of f\ on < 0.2 and f mn > 1.0. The peak at xm = 1 is over 3.4 times as 
likely as the next closest value of xm = 0.17, in which 7% of the LOSs pass. 

3. RESULTS 

The procedure outlined above turns out to provide tight con- 
strains on all three of our free parameters simultaneously. The 
results are summarized in Figure 2, showing the percentage of 
LOSs that pass our analysis criteria, as a function of Rs, xui, 
and /{ 0n . The fraction of accepted LOSs peaks sharply at a 
maximum of 24% at the parameter choice of (Rs, xhi, /ion) 
= (42 Mpc, 1, 0.4). The radius of the ionized sphere is lim- 
ited to 42 Mpc < Rs < 47 Mpc, with no LOSs passing our 
criteria outside of that range. Note that an offset in the Ly- 
man a line center by ±1000 km s" 1 (typical of quasar jets, 
for example), would represent an additional ±1.5 Mpc uncer- 
tainty on the inferred radius. Although SDSS J1030+0524 has 
a precise redshift determination from C IV and N V lines, a 



redshift error of Az = 0.01-0.02 (typical in the absence of 
such metal lines) would add an uncertainty of ±(4-8) Mpc 
to any analysis similar to the one presented here. The ion- 
ization flux parameter is tightly constrained around fi on = 0.4, 
with no LOSs passing our criteria with values of fi on < 0.2 and 
/i on > 1 .0. The peak value of / on = 0.4 implies an emission rate 
of (5.2 ± 2.5) x 10 56 photons sec" 1 . Finally, the peak at xhi = 1 
is over 3.4 times as likely as the next closest value of xm = 0.17, 
in which 7% of the LOSs pass, and 6 times as likely as values 
of xhi < 0.016, in which 4% of the LOSs pass. 

These results can be interpreted as follows. As mentioned 
previously, the presence of flux in Region 2 (tl yq < 22.8) sets a 
lower limit on Rs, implying that the edge of the ionized region 
must be at a smaller wavelength than the boundary between Re- 
gions 2 and 3. Region 3, however, yields an upper limit on Rs, 
from the requirement that TL ya > 22.8 in that region. This high 
optical depth cannot be maintained by tr alone, without vio- 
lating the constraint in Region 1 of T\^ a < 6.3. We emphasize 
that the upper limit on Rs makes no use of the data in Region 
2. Hence the edge of the ionized region must be close to the 
boundary between Regions 2 and 3. The ionization flux is con- 
strained due to our criteria in Regions 1 and 2. An ionization 
flux that is too small (large tr) would violate TL ya < 6.3 in Re- 
gion 1 . An ionization flux that is too large (small tr) would vi- 
olate TLy Q > 6.3 throughout Region 2, since the damping wing 
contribution to tl yq is insufficient, given the lower limit on Rs 
above. Our constraint on the neutral hydrogen fraction, xm, 
comes from the presence of flux in the Lyman j3 region of the 
spectrum corresponding to Region 2. Because of fluctuations 
in the density field (and hence in t«), a strong damping wing is 
needed to raise 7L yQ above 6.3 throughout Region 2, while still 
preserving Ti^ a < 6.3 in Region 1. 

4. DISCUSSION AND CONCLUSIONS 

Using only the lower limits on the Lyman a optical depth 
obtained from the Lyman a and Lyman f3 GP troughs of SDSS 
J 1030+0524, we find a simultaneous constraint on the size of 
the ionized region surrounding this quasar, the neutral fraction 
of the ambient intergalactic medium, and the ionizing luminos- 
ity of the quasar. 

The radius of the Stromgren sphere is limited to 42 Mpc 
< Rs < 47 Mpc, close to the previously estimated lower lim- 
its (Madau & Rees 2000; Cen & Haiman 2000). We also find 
evidence that the IGM was significantly neutral at z ~ 6, with a 
~ 1 a lower limit of xhi > 0.17. This result is derived from the 
observed sharpness of the boundary of the HII region alone, and 
relies only on the gross density fluctuation statistics from the 
numerical simulation. In particular, it does not rely on assump- 
tions about the mechanism for the growth of the HII region. 
Our results provide a robust confirmation of previous sugges- 
tions that the IGM was significantly more neutral at z ~ 6 than 
the lower limit of xm > 10~ 3 that is directly obtainable from the 
black region (Gunn-Peterson trough) of the quasar spectra. Fi- 
nally, we find an emission rate of ionizing photons per second 
of (5.2 ± 2.5) x 10 56 photons sec" 1 , which is between 2-10 
times lower than expected (Elvis et al. 1994; Telfer et al. 2002), 
strengthening arguments that reionization at z ~ 6 is caused by 
the radiation from early stars, rather than bright quasars (e.g. 
Dijkstra et al. 2004). 

Our findings represent the first detection of the boundary of 
a cosmological HII region, and have several important impli- 
cations. It has recently been shown that tight constraints on 



the hydrogen neutral fraction can be extracted directly from the 
Lyman a absorption spectrum alone (Mesinger et al. 2004), but 
this requires tens of spectra to be statistically significant when 
Rs is not known. Incorporating an independent limit on Rs, 
such as the one obtained from this method, can reduce the num- 
ber of required spectra to one. Tight constraints on xhi have also 
recently been obtained by adopting a proper size of Rs ~ 4.5 
Mpc for SDSS J1030+0524, together with an estimate of its 
lifetime (Wyithe & Loeb 2004). Our direct determination of the 
Stromgren sphere size is only slightly larger than the previously 
adopted value, lending further credibility to this conclusion. 

The sharp boundary we detect also constrains the hardness 
of the ionizing spectrum of SDSS J1030+0524. We infer here 
a rise in the neutral fraction over a redshift interval of Az ~ 
0.02, corresponding to a proper distance of ~ 1.2 Mpc. In 
order for its mean free path not to significantly exceed this 
thickness, the energy of the typical ionizing photon emitted by 
SDSS J1030+0524 must be < 230eV. This implies that the ef- 
fective slope of the ionizing spectrum of this source is softer 
than-<f lnL^/if lnz/= 1.07. This analysis ignores radiative trans- 
fer effects, which can further blur the apparent boundary of the 
HII region, and would strengthen this limit. 

The analysis presented here uses only the Lyman (3 and Ly- 
man a absorption spectra for the single quasar J 1030+0524, 
but can be extended by incorporating additional Lyman lines, 
which provide further constraints on the optical depth. Al- 
though such analysis will have additional uncertainties associ- 
ated with the need to subtract foreground absorption due to the 
lower Lyman transitions, it is likely that adding the Lyman 7 
line can tighten the constraints obtained here (Fan et al. 2002). 
Note that the presence of dust within the HII region can atten- 
uate and redden the observed spectrum. This will not, how- 
ever, impact our conclusions, since the overall attenuation is 
absorbed into /j on , and reddening extends over a wavelength 
range much broader than the sharp features considered here. 

In the future, given a larger sample of quasars (and/or a sam- 
ple of gamma-ray burst afterglows with near-IR spectra) at z > 
6, it will be possible to use the method presented here to search 
for sharp features in the absorption spectrum from intervening 
HII regions, not associated with the background source itself. 
We plan study the utility of this approach in a future paper. For 
such external HII regions that happen to intersect the line of 
sight, both the red and the blue sides of each GP trough can, 
in principle, be detected, and can yield two separate measure- 
ments of the ionized fraction at different points along the line of 
sight. The Universe must have gone through a transition epoch 
when HII regions, driven into the IGM by quasars and galax- 
ies, partially percolated and filled a significant fraction of the 
volume. The detection of the associated sharp features in fu- 
ture quasar absorption spectra will provide a direct probe of the 
the 3D topology of ionized regions during this crucial transi- 
tion epoch. Deep surveys equipped with sufficiently red filters 
(with instruments such as those being carried out by the VLT, 
and ultimately with deep surveys covering a significant portion 
of the sky, such as the survey proposed with the Large-aperture 
Synoptic Survey Telescope, LSST; Tyson (2002)) will be able 
to deliver a large sample of bright z > 6 quasars needed for such 
studies. 
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